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The approach is developed for the description of isolated Fermi-systems with finite number of particles, such as complex 
atoms, nuclei, atomic clusters etc. It is based on statistical properties of chaotic excited states which are formed by the 
interaction between particles. New type of "microcanonical" partition function is introduced and expressed in terms of the 
average shape of eigenstates F{Ek,E) where E is the total energy of the system. This partition function plays the same 
role as the canonical expression ea;p(— i?'''/r) for open systems in thermal bath. The approach allows to calculate mean 
values and non-diagonal matrix elements of different operators. In particular, the following problems have been considered: 
distribution of occupation numbers and its relevance to the canonical and Fermi-Dirac distributions; criteria of equilibrium and 
thermalization; thermodynamical equation of state and the meaning of temperature, entropy and heat capacity, increase of 
effective temperature due to the interaction. The problems of spreading widths and shape of the eigenstates are also studied. 
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I. INTRODUCTION 

As is known, quantum statistical laws are derived for 
systems with infinite number of particles, or for systems 
in a heat bath, therefore, their applicability to isolated 
finite systems of a few particles is, at least, question- 
able. However, the density of many-particle energy levels 
increases extremely fast (typically, exponentially) both 
with an increase of number of particles and excitation 
energy. For this reason, even a weak interaction be- 
tween particles can lead to a strong mixing between large 
number of simple many-particle states, resulting in the 
so-called chaotic eigenstates. If the components of such 
eigenstates can be treated as random variables (onset of 
Quantum Chaos ), statistical methods are expected to 
be valid even for an isolated dynamical system. 

One should stress that statistical description of such 
isolated systems can be quite different from that based 
on standard canonical distributions; therefore, applica- 
tion of the famous Fermi-Dirac or Bose-Einstein formu- 
lae may give incorrect results. Moreover, for isolated 
few-particle systems a serious problem arises in the defi- 
nition of temperature, or other thermodynamic variables 
like entropy and specific heat (to compare with, in infi- 
nite systems different definitions give the same result). 

The aim of this paper is to develop statistical theory for 
finite quantum systems of interacting particles, based on 
generic statistical properties of chaotic eigenstates ("mi- 
crocanonical" approach). Typical examples of such sys- 
tems are compound nuclei, complex atoms, atomic clus- 
ters, isolated quantum dots, etc. 

The structure of the paper is as follows. In the next 
Section 2, new type of "microcanonical" partition func- 
tion is introduced for finite isolated systems, which is 
directly related to the average shape of chaotic eigen- 
states. Based on this partition function, general expres- 
sion for the occupation numbers is given which is valid 
for any number of interacting particles. Relation of this 



"microcanonical" expression {F — distribution) to that of 
standard canonical distribution for occupation numbers 
is the subject of the Section 3. Also, here a new form 
of the canonical distribution is suggested which can be 
convenient in applications to open systems in thermal 
equilibrium like quantum dots. 

Transition to the standard Fermi-Dirac distribution for 
large systems is analyzed in Section 4. Specific attention 
is paid to the accuracy of this distribution in applications 
to isolated systems with few particles. 

In Section 5 the analytical derivation of the 
-F— distribution for occupation numbers is given. For 
this, the model of n randomly interacting Fermi-particles 
distributed over m single-particle levels has been used. 
The analytical formula for the occupation numbers Us (E) 
with E being the total energy of the system, was found to 
be in excellent agreement with the numerical experiment. 

In Section 6 the influence of finite width of single- 
particles states ( "quasi-particles" ) on the distribution of 
occupation numbers is considered. New expression for 
the occupation numbers is discussed which takes into ac- 
count the finite spreading width of "quasi-particles" . It is 
demonstrated that for an isolated system with the fixed 
total energy E , the incorporation of spreading widths 
decreases the effective temperature, however, does not 
change the occupation numbers. 

Section 7 deals with an important problem of ther- 
modynamical description of small systems consisting of 
finite number of interacting particles. Specifically, dif- 
ferent definitions of both temperature and entropy have 
been analyzed, and the equation of state for finite sys- 
tems has been derived. 

In Section 8 we show that statistical effects of the in- 
teraction can be imitated by an increase of the effective 
temperature. This fact allows to use the standard Fermi- 
Dirac expression for the occupation numbers (with renor- 
malized parameters) in the application to both isolated 
and open (in the thermal bath) systems of interacting 
Fermi-particles. This interaction may be strong, how- 
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ever, the assumption of randomness for residual interac- 
tion matrix elements is essential. 

In Section 9 the conditions for chaos, equilibrium 
and "thermalization" have been analyzed for an isolated 
quantum system of a finite number of particles in terms 
of the interaction strength, number of particles and other 
related parameters. Depending on these conditions, there 
are four different regions with different distributions of 
occupation numbers, which are discussed in detail. In 
this Section, the transition to mesoscopic systems is 
shortly discussed. The question of particular interest is 
how statistical properties depend on the dimensionality 
of a system. 

Since the approach developed in this paper is entirely 
related to the structure of chaotic eigenstates, in Ap- 
pendix 1 the analytical expression for the average shape 
of eigenstates is given and discussed in more details. This 
expression, which also describes the local spectral density 
of states (LDOS), essentially depends on two different 
widths and is valid both for weak and strong interac- 
tion between particles. For weak interaction the shape is 
close to the Lorentzian form with the half- width given by 
the Fermi Golden Rule. For larger interaction the shape 
is characterized by exponential tails and by the width 
which is linear in the interaction strength (contrary to the 
Fermi Golden Rule which gives a quadratic dependence). 
Transition between these two regimes occurs when the 
half-width is comparable to the root-mean-square width 
(effective band width of the Hamiltonian matrix). 

In Appendix 2 the moments of the distributions of the 
basis components (LDOS) and energy levels are calcu- 
lated. 

And finally, in Appendix 3 the calculation of density of 
final states and spreading widths of the LDOS has been 
performed using the Fermi Golden Rule. 



II. MICROCANONICAL PARTITION FUNCTION 

In this Section we are going to derive the partition 
function for closed (isolated) systems of finite number of 
interacting particles. This function allows one to per- 
form analytical and numerical calculations of statistical 
mean values of different operators, for example, occupa- 
tion numbers. 

We follow the standard approach which is based on the 
separation of a total Hamiltonian in two parts, 

H ^ Ho + V e^a+fls + ^ X! ^pgrsOp a^flrfls (1) 

The "unperturbed" Hamiltonian Hq should incorporate 
the effect of a mean field (if any), are the energies of 
single-particle states ( "orbitals" ) calculated in this field, 
aj",as are creation-annihilation operators, and V stands 
for the residual interaction. For simplicity, here we ne- 
glect any dynamical effects of the interaction like pairing, 
collective modes etc. Instead, we will study statistical ef- 
fects of interaction, therefore, in what follows we assume 



that matrix elements Vpqrs of the two-body residual in- 
teraction are random variables. 

Exact ("compound") eigenstates \i) of the Hamiltonian 
H can be expressed in terms of simple "shell-model basis 
states" I fc) (eigenstates of i/o) : 
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\k); \k)^. 



(2) 



These compound eigenstates \i) , characterized by the 
corresponding energies i?*-'-* , are formed by the residual 
interaction V ; in complex systems they typically contain 
large number Npc ^ 1 of the so-called principal compo- 
nents C^*' which fiuctuate "randomly" as a function of 
indices i and k . 

Our main interest is in the occupation numbers of 
single-particle states (orbitals). They can be represented 
in terms of components of the exact eigenstates \i) , 



(fc| fis \k) 



(3) 



Here fis — a^ag is the occupation number operator. 
The knowledge of the distribution of occupation num- 
bers gives the possibility to calculate mean values of any 
single particle operator (M) = "^sMss ■ Moreover, 
the variance of the distribution of non-diagonal elements 
of M describing transition amplitudes between "chaotic" 
compound states, can be also expressed through the oc- 
cupation numbers . 

As one can see from (3), mean values of occupation 
numbers depend on the shape of exact eigenstates, given 
by the "spreading function" F (in what follows, the F — 
Junction ), 



F, 



F(i?fc-i?«) 



(4) 



Last equality in the above expression reflects the fact that 
residual interaction V mainly mixes close components k 
in some energy interval F ( "spreading width" ) nearby 
the eigenstate energy E^"^^ (more accurately, nearby the 
unperturbed energy Ek for k = i , see below) . 

Typically, this spreading function rapidly decreases 
with an increase of | — E^^^^ | (since an admixture of 
distant component is very small) . Recently, in numerical 
studies of the Ce atom [g] , the s — d nuclear shell model 
and random two-body interaction model |3|j|] it was 
found that typical shape F of exact eigenstates practi- 
cally does not depend on a particular many-body system 
and has a universal form which essentially depends on 
the spreading width F . The latter can be expressed in 
terms of parameters of the model ( intensity V of the 
residual interaction, number of particles n , excitation 
energy, etc.) and can be calculated analytically (see Ap- 
pendices 1- 3). One can also measure the width of the 
_F— function (0) via the number of principal components, 
Npc ''^ T/D where D is the local mean energy spacing 



2 



for compound states. In many-body systems the value of 
D exponentially decreases with an increase of the num- 
ber of "active" (valence) particles. As a result, Npc is 
very large, ~ 10^-;- 10^ in excited (compound) nuclei and 
~ 100 in excited rare-earth or actinide atoms. 

The starting point of our approach is the expression for 
the occupation numbers which stems from Eqs. (^,|^), 



n,{E) 



Y^nf^ F{Ek~E) 

_k 

Y.F{Ek-E) 

k 



(5) 



where rS^^ = {k\hs\k) equal or 1 for Fermi- 
particles and the sum in the denominator stands for the 
normalization. This way of averaging of occupation num- 
bers is a kind of microcanonical averaging since it is de- 
fined for the fixed total energy _E of a system. One can 
see that, in fact, the relation (j|) is equivalent to the in- 
troduction of a new kind of partition function. 



(6) 



which is entirely determined by the shape of chaotic 
eigenfunctions. In what follows, we term Eq.(^) the 
F — distribution. 

The above expression (|]) gives a new insight on the 
problem of statistical description of complex systems. In- 
deed, as was mentioned above, the shape of the F— func- 
tion has universal features and can be often described an- 
alytically, therefore, in practice there is no need to diago- 
nalize huge Hamiltonian matrix of many-body system in 
order to find statistical averages. One should stress that 
the summation in (^ is carried out over unperturbed en- 
ergies Ek defined by the mean field, rather than over the 
energies of exact eigenstates in the standard canonical 
distribution. As a result, the distribution of occupation 
numbers can be derived analytically (see Section 5) even 
for few interacting particles, in the situation when the 
standard Fcrmi-Dirac distribution is not valid. 



III. TRANSITION TO THE CANONICAL 
DISTRIBUTION 

It is instructive to compare our F— distribution (|^) 
with occupation numbers obtained by making use of the 
standard canonical distribution, 



ns{T) 



Y.n'i^ exp(-i;(*)/T) 

i 



(7) 



where T is the temperature and the index i stands for 
exact eigenstates. The important difference between the 
F— distribution (||) and the canonical distribution (|^) is 
that in Eq. (|^) the occupation numbers are calculated 



for a specific energy i? of a system unlike specific tem- 
perature T in Eq.(0). However, results of calculations 
based on Eqs. (^ and (||) can be compared with each 
other using the relation between the energy E and the 
temperature T , 



E = (E), 



cxp(-£;W/T) 

i 



(8) 



The comparison of Eqs. ^ and (||) also shows that the 
canonical distribution corresponds to the averaging of the 
"microcanonical" F— distribution over some energy in- 



terval At . To demonstrate this, let us substitute n, 
and C, 



(0 



from Eqs.(^,^ into Eq.(|^) and replace the 

summation over i by the integration over p{E''^^)dE'^'^^ 
where p^E^"^^ ) is the density of exact energy levels. 



/"nW$2,(i;W)d£;W (9) 



n^'^ exp I 



Here we have introduced the "canonical (thermal) aver- 
aging" function, 



^t{E) = p(£;)exp i~E/T) 



(10) 



which is discussed below. As a result, we can transform 
canonical distribution (|^) into the form similar to the F~ 
distribution (O), 



risiT) 



k 



(k) 



F{T,Ek) 



T.F{T,Ek) 

k 



(11) 



where the function F{T,Ek) is the canonical average of 



F, 



(0 



F{T,Ek) - J Fi^^^T{E^-^'>)dE'^^ 



(12) 



Note that this form of the canonical distribution can 
be convenient for the calculation of the occupation num- 
bers and other mean values in quantum dots which are in 
thermal equilibrium with an environment (with no par- 
ticle exchange). 

In large many-body systems the canonical averaging 
function (^t{E) has narrow maximum since the density 
of states p{E^^'>) typically grows very fast. The position 
Em of its maximum is defined by the expression 



d\np{E) _ 1 
dE ^ T 



and the width is given by 



^2 hip{E) 



dE^ 



-1/2 



(13) 



(14) 
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As an example, let us consider the system of n interacting 
particles distributed over m orbitals. In the papers [|6|J3] 
it was shown that in the case to 3> n 3> 1 the density of 
states is of the Gaussian form, 



p{E) = 



1 



exp 



2a2 



(15) 



where Ec is the center of the spectrum and is its vari- 
ance. According to recent numerical data , the Gaus- 
sian form for the density p{E) occurs also for few parti- 
cles (n > 4) . This fact allows easily to find the form of 
^t{E) which appears to be of quite generic. 



$t(-B) ^ exp 



{E ~ E^f 
2a2 



where 



Er, 



Ec 



T 



(16) 



(17) 



One can see that the width Ar of the thermal averaging 
function equals to the root-mean-square (r.m.s.) width 
of the spectrum, A-p = a . Now, it is easy to show that 
the thermal averaging width Ay is always bigger than the 
r.m.s. width AE' of the "microcanonical" F— function or, 
the same, than the mean width of exact eigenstates in the 
energy representation. Indeed, there is a simple relation 
between the widths a and (Tq of the energy spectrum with 
and without interaction, respectively, (see Appendix 1), 

= al + {^Ef , (18) 

therefore, we have At = (J > AE . One should stress 
that the latter width AE , in fact, is due to statistical 
effects of interaction. The difference between the widths 
At and AE is not important when the number of par- 
ticles n is large. This is because with an increase of n 
the width ctq of the unperturbed spectrum increases as 
•y/n, unlike the width AE which increases as n. One 
should also note that in this case both widths At and 
AE are much smaller than the typical energy interval, 
a/\E — Ec\ ~ • Therefore, for large number of 

particles the function $t can be regarded as the delta- 
function ai E = Em and the F— distribution is close to 
the canonical distribution, see Eq. (p^ . 

To conclude with this Section, the canonical distribu- 
tion (^ is not correct when describing isolated systems 
with small number of particles, instead, one should use 
the F— distribution (|g). This was recently confirmed 
by numerical experiments with the model of few Fermi- 
particles with two-body random interaction MM. 



IV. TRANSITION TO THE FERMI-DIRAC 
DISTRIBUTION 

It is instructive now to show how the standard Fermi- 
Dirac distribution stems directly from the F— distribu- 
tion (ph in the limit of large number of particles. By 



performing explicitly the summation over = and 
Us = 1 , the expression (^ can be rewritten in the form 



n,{E) 



+ Z,{n-l,E-es) ^ _ 
Z,{n-l,E-is) + Z,{n,E) i 



1 



Z,{n,E) 

(19) 



Here, two "partial" partition functions Zs{n,E) and 
Zs{n — 1,E ~ is) are introduced. In the first one, the 
summation is taken over all single-particle states of n par- 
ticles with the orbital s excluded, Zs{n, E) — ^'j. F{Ek — 
E). Correspondingly, the sum in Zg{n—\, E—is) is taken 
over the states of n — 1 particles with the orbital s ex- 
cluded. The latter sum appears from the terms for which 
the orbital s is filled {ug = 1) , thus, we should add the 
energy is = Ek{n) — Ek{n — 1) of this orbital to the en- 
ergy Ek{n — 1) of the basis state \k) defined by n — 1 
particles. Since the F— function depends on the differ- 
ence Ek + is — E only, the adding term to Ek{n — 1) 
is the same as its subtraction from the total energy E . 
Note, that this term is defined by 



(20) 



where ts is the energy of a single-particle state and Ugp 
is the diagonal matrix element of the two-body residual 
interaction. By taking independent of k we assume 
that the averaging over the basis states near the energy 
E is possible, in fact, this is equivalent to a local (at a 
given energy) mean field approximation. 

One should stress that this approximation is the most 
important when applying the model (0) to realistic sys- 
tems. For example, for Ce atom there are orbitals from 
different open sub-shells (e.g. 4/ and 6s ) which are 
quite close in energies, however, they have very different 
radius. As a result, the Coulomb interaction between the 
corresponding electrons is very different In this case 
the interaction terms in Eq.(|2C|) strongly depend on the 
occupation numbers of other particles, which means that 
there is no good mean field approximation. As a result, 
the equilibrium distribution for occupation numbers is 
very different from the Fermi-Dirac distribution |^ . How- 
ever, the F— distribution (^) for occupation numbers is 
valid. In other cases like random two-body interaction 
model or nuclear shell model [|| , such a local mean 

field approximation is quite accurate. 

For large number n 3> 1 of particles distributed over 
TO ^ 1 orbitals, the dependence of Zg on n and is is 
very strong since the number of terms N in the partition 



function Zs is exponentially large, N 



There- 



(m — n) !n! 

fore, to make the dependence on arguments smooth, one 
should consider In instead of Zg ■ In this case one can 
obtain 



In Zsin — An, E ~ ig) ^ In Zgin, E) — oig An — (3gig 

An = 1 (21) 



dlnZg ^ dlnZg 
ois = — ; Ps = 



dn 



dE 
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This leads to the distribution of the Fermi-Dirac type, 



1 



1 + exp(as + Pais) 



(22) 



If the number of substantiaUy occupied orbitals in the 
definition of Zs is large, the parameters as and fig are 
not sensitive as to which particular orbital s is excluded 
from the sum and one can assume as = a = — /z/T, (3s = 
(3 = 1/T as in the standard Fermi-Dirac distribution. 
Then, the chemical potential fi and temperature T can 
be found from the conditions of fixed number of particles 
and fixed energy, 



J2 ^s-ris + J2 



Y.ns[es + es)/2^E 



(23) 



Note, that the sums in (^ , ^0|) containing the residual 
interaction Usp can be substantially reduced by a proper 
choice of the mean field basis (for instance, the terms Usp 
can have different signs in such a basis). In practice, the 
values ts and Eg may be very close. Since in the above 



expressions (23) the non-diagonal matrix elements of the 
interaction are not taken into account, one can expect 
that the distribution of occupation numbers defined by 
these equations gives a correct result if the interaction 
is weak enough (ideal gas approximation). However, we 
can show (see Section 8) that, in fact, even for strong 
interaction the Fermi-Dirac distribution can be also valid 
if the total energy E is corrected in a proper way, by 
taking into account the increase of the temperature due 
to statistical effects of interaction. 

One should also note that somewhat similar procedure 
transforms the canonical distribution (^ into the Fermi- 
Dirac distribution (see e.g. |l^) in the case of many non- 
interacting particles (ideal gas). It is curious that the 
Fermi-^Dirac distribution is very close to the canonical 
distribution (^ even for very small number of particles 
{n = 2), provided the number of essentially occupied 
orbitals is large (which happens for T 3> e or /i 3> e ) . In 
fact, this is a result of a large number of "principal" terms 
in the partition function Zs which allows us to replace Us 
by a in the term Zs{n,T)/Zs(n — 1,T) = exp(as -I- /3T) 
in the canonical distribution (^ (compare with (p9|)). 

More accurate consideration shows, however, that the 
temperature T in the Fermi-Dirac distribution is differ- 
ent from that in the canonical distribution. Indeed, using 
the expansion as — alep) + a'{es — ep) one can obtain 
the relation between the Fermi-Dirac (Pfd) and canoni- 
cal (P) inverse temperatures, Pfd — fS + a'ep ■ Concern- 
ing the chemical potential, its definition also changes, 
~fJ./T — alep) — a'ep . This fact is confirmed by our 
numerical simulations for an isolated system with few in- 
teracting Fermi-particles ||,^,||. Namely, for the same 
total energy E of the system , the canonical and Fermi- 
Dirac distributions give the same distribution Us defined, 
however, by different temperatures since they are deter- 
mined by different equations (||) and (p3|). 



The closeness of these two distributions for any num- 
ber of particles is not so surprising in the presence of the 
thermostat, where even one particle is in the equilibrium. 
To the contrary, for isolated systems with small number 
of particles the applicability of the Fermi-Dirac distribu- 
tion is not obvious. To answer this question, one needs 
to analyze the role of interaction in the creation of an 
equilibrium distribution. 



V. ANALYTICAL CALCULATION OF 
OCCUPATION NUMBERS IN FINITE SYSTEMS 

The advantage of the approach developed in this pa- 
per is that if we know the shape of eigenstates in the 
many-particle basis of non-interacting particles (the F— 
function) and the unperturbed density of states po{E), 
one can analytically calculate the distribution of the oc- 
cupation numbers Ug. 

In order to calculate the occupation numbers n^, we 
use the expression ( p^ ) containing two partial partition 
functions Zs{n, E) and Zs{n—1, E—es) which correspond 
to systems with n and n — 1 particles, with the orbital 
s is excluded from the set of single-particle states. The 
partition function can be found from the relation 



Z = Y,FiEk~E)^ f po{Ek)F{Ek - E)dEk 

I, 



(24) 



The density of unperturbed states po{Ek) in a system 
of n particles distributed over m single-particle states 
(orbitals) was shown to be close to the Gaussian (see for 
example, [|i,lll): 



N 



exp 



{Ek — Ec 



(25) 



with Ec as the center of the energy spectrum and N as 
the total number of states. Let us assume that the shape 
of eigenstates F— is given by the Gaussian, too: 



F{Ek -E) = 



1 



exp 



{Eu - Ef 
2(A£;)2 



Here the variance (Ai?)^ is defined by Eqs. 



(26) 



and 



E = E^^ + A 



(i) 



Ha (see Appendix 1). 



As was found, Gaussian shape of eigenstates (apart 
from long tails) occurs in realistic systems like Ce atom 
1^ and heavy nuclei . Recently, the form of the 
F— function in dependence on the perturbation has been 
studied in details in the model of Wigner Band Ran- 
dom Matrices , as well as in the random two-body inter- 
action model [|| (see also nuclear calculations p^ , p^ ). 
In particular, it was discovered, that the gaussian-type 
shape happens when the interaction is large enough, 
namely, when the Breit- Wigner width F — 2i:pV^ is com- 
parable with the root-mean-square width Ai? (an effec- 
tive band width of a Hamiltonian matrix) . 
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By performing the integration in Eq.(p4[) one gets 

\2 \ 



Z[E) 



N 



exp 



(27) 



where cr^ = ctq + {/S.EY (it coincides with the variance of 
the perturbed spectrum). In order to calculate the occu- 
pation numbers ris, we use the expression (19). For this, 
one needs to find the partial partition functions Zs(n, E) 
and Zs{n — 1, E — es) corresponding to n and n~l parti- 
cles with the orbital s is excluded from the set of single- 
particle states. Now we have to calculate the number of 
states Ng and the center Ecs for these truncated systems, 



iV,(n,m- 1) 



Ns{n~ l,m- 1) 



(m- 1)! 



(to — 1 — n)! nl 



(to - 1)! 



(m — n)! (ri — 1)! 



Ecs{n) = e-s n ; £'cs(n - 1) = {e-s){n - 1) 



_ E 



TO — 1 



The variance ctos of the energy distributions can be esti- 
mated as 



where afg is the variance of single-particle spectrum with 
the excluded orbital s . Here, for simplicity, we have 
neglected the Pauli principle which is valid for m ^ n . 
More accurate calculation can be easily done with the use 
of calculator. As a result, the distribution of occupation 
numbers has the form 



rniE) 



1 



1 + R 



R 



m — n as{n — 1) 
n as (n) 



X exp 



iE-E,s{n)f ^ [E - es ~ E,s{n - l)f 



2al(n) 



2a2(n-l) 



(28) 



where al = cr^g + {^Ef . 

It is instructive to compare this result with the Fermi- 
Dirac distribution which is valid for large number of 
particles. In this case R = exp((es — ij)/Tth) where 
Tth = cf'^/{Ec — E) is the thermodynamic temperature 
which is discussed below, see (|33|). The chemical poten- 
tial jj, in this case should be calculated numerically to fix 
the total number of particles n. The data are reported 



in Fig.l. One can see that Eq.(^8|) predicts occupation 
numbers in perfect agreement with the numerical exper- 
iment. 

Finally, note that the same method can be used to 
solve another problem: to find the distribution of the oc- 
cupation numbers ns{T) in finite systems of interacting 
particles in the thermal bath with the temperature T . 
For such a case, it is enough to replace the i^— function 
by the canonical average F{T,Ek), see (|l|). In fact, it 
is the method for taking into account the "random" in- 
teraction in the canonical distribution. The result for 
the occupation numbers Ug (T) can be obtained from ( [2^ ) 
by replacing a"^ 2cr^ and E — Em + Ai{Em) where 
Em ^ Ec — a"^ /T and Ai is a small correction, see Ap- 
pendix 1. 



VI. PARTICLES AND QUASI- PARTICLES, ROLE 
OF SINGLE-PARTICLE SPREADING WIDTH 

In previous sections we have discussed the distribu- 
tion of occupation numbers for real particles distributed 
over given orbitals. At the same time, there exist tradi- 
tional approach which is based on the notion of "quasi- 
particles" . It allows to incorporate the effects of interac- 
tion in terms of single-particle states and goes beyond the 
mean-field approximation. As is well known, the inter- 
action leads to the spreading width 7s for single-particle 
orbitals . It also results in the shift of average energies, 
Es = Es -I- 5ts ■ According to our numerical data ( j|]) for 
the random two-body interaction model, the shifts ^e^ 
turn out to be, in average, smaller than 7^ and for this 
reason one can take into account the effect of spreading 
widths 7s only. 

Here we would like to analyze the role of the spreading 
widths for the distribution of occupation numbers and 
compare with our approach where the interaction is taken 
into account in terms of many-body states. For this, let 
us average the standard Fermi-Dirac occupation numbers 
Us over the energy interval 7s: 



-Is/2 



n{e) — = 1 In 

Is 7s 



1 -I- exp 



2T 



_ 1 -I- exp ■ 



2T 



where 



n(e) = 



l + cxp(^) 



(29) 



(30) 



It seems, this is the simplest form of the Fermi-Dirac 
distribution for the "quasi-particles" with finite spread- 
ing widths. One can check that in the limit 7s = the 
Fermi-Dirac expression ( pO[ ) with Ug = n{es) is recovered. 

To test the sensitivity of the occupation numbers to the 
values of the spreading widths 7s, we have solved equa- 
tions (23) for chemical potential and temperature (for a 
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given energy of our isolated system) using the standard 
expression (^0|), and compared the result with that ob- 
tained by using the expression (p9|), see details in 
The data have revealed that the chemical potential prac- 
tically does not change, while there is a noticeable de- 
crease of the temperature, T(7 7^ 0) < T{j = 0) . The 
striking result is that the two curves for the occupation 
numbers (^9|) and ( |30| ) coincide with a high accuracy, 
namely, ns(es,7,T) sa 713(6,5,7 = 0,T + AT) = ns{es,f). 
This means that the temperature mimics the effect of the 
spreading widths, the phenomenon which is far from be- 
ing trivial. The shift of the temperature for 7 <C /i can 
be estimated analytically as AT « 7^/16T. 

The above result indicates that one should not worry 
about the finite width of single-particle orbitals ( "quasi- 
particles") when calculating the occupation numbers. 



VII. THERMODYNAMICS OF SMALL SYSTEMS 

One of the important questions is about thermody- 
namical description of isolated systems of interacting par- 
ticles. In any thermodynamical approach one needs to 
define, in a consistent way, such quantities as entropy, 
temperature and equation of state. Different definitions 
of the entropy and temperature have been recently dis- 
cussed ll^] (see also [0) in application to shell models 
of heavy nuclei. In particular, it was found that for a 
realistic residual interaction, different definitions of tem- 
peratures lead to the same result. Below, we analyze few 
definitions of the temperature and entropy and show that 
for small number of interacting particles they may give 
quite different results. 

Standard thermodynamical definitions of the entropy 
Sth and temperature are based on the density of states 



Sth = In p{E) + const 



1 

Tth 



dSth d\n p 



dE 



dE 



(31) 



(32) 



In fact, such a definition of the temperature follows from 
the estimate of the position of maximum of the canonical 
averaging function <i>T(£') , see ( p^ and (p^). It is usually 
assumed that the position of its maximum Em coincides 
with the energy £^ of a system. One should stress that 
in the above definitions p{E) is the perturbed density of 
states, therefore, the interaction is essentially taken into 
account. 

However, for finite isolated systems with the fixed en- 
ergy E , the definition of the temperature given by the 
relation {E)rp = E (see Eq.(H) seems to be more natu- 
ral. Here the averaging is performed over the canonical 
distribution (Q). Since the width At of the canonical av- 
eraging function '^t(E) is not zero, the two definitions of 
the temperature, ( ^2|) and (|^) give, in principal, different 



results. Indeed, in the case of the Gaussian form of p{E) 
the value of Tth given by (^2|) takes the form (see also 

0), 



Tth — 



Ec-E 



(33) 



where E^ and a are the center and the width of the dis- 
tribution p{E). 

On the other hand, direct evaluation of the relation (||) 
leads to the following definition of the temperature. 



(34) 



Here, the shift A is given by the expression 



(35) 



where 



K = 



exp 



-^\dx 



/2n 



E — Er, 



E„ 



Er 



Tr. 



(36) 



One can see that the shift A itself depends on the tem- 
perature and is proportional to the width At =<j of the 
function $t(^')- In the above relations, i?min and -Emax 
are the low and upper borders of the energy spectrum. 
Note that the relation A = occurs at the center of 
the spectrum, therefore, the temperature in the upper 
part of the spectrum is negative (it is typical for systems 
with bomided spectrum, for example, for spin systems). 
In fact, our model (^ with finite number ra of orbitals 
can be treated as the model of one open shell in atoms, 
nuclei, clusters, etc. However, in realistic many-body sys- 
tems there are always higher shells which contribute to 
the density of states for higher energy. Thus, the density 
of states p{E) is a monotonic function which results in 
the positive temperature. For such physical applications, 
the model ([I]) with finite number of orbitals is reasonable 
in the lower part of the energy spectrum where the influ- 
ence of higher shells can be neglected. 

One can also see that the difference between the 
two equations of state T{E) defined by Eqs.(^) and 
( p^ ) , disappears for highly excited eigenstates (for which 
— T^min » cr), or in large systems with 1 . Indeed, 
one can obtain, E^ — E nai , where ai is the width of 
single-particle spectrum. On the other hand, according 
to the central limit theorem, the variance of total energy 
spectrum can be estimated as ctq « X^n ""i = ''^'^1 ' there- 
fore, the ratio a/{Ec—E) ~ l/v^ tends to zero at n — > 00 
. Note, that in finite systems (atom, nucleus etc.) the 
number of valence particles (particles in an open shell) is 
not large. For example, for Ce atom we have rt = 4 ||] 
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and in nuclear shell model [Q n = 12 , therefore, the cor- 
rections to the thermodynamical temperature ( ^2|) can 
be significant, especially, for low energies. Here, we do 
not take into account particles from deep closed shells 
since their excitation energy is high and they do not con- 
tribute to the thermodynamical and statistical proper- 
ties of systems (though they renormalize parameters of 
the Hamiltonian (|l|) describing the interaction between 
valence particles). 

The energy dependence of temperatures Tth and Tcan is 
shown in Fig. 2. The data are given for the model oin = 4 
interacting Fermi-particles distributed over m = 11 or- 
bitals. The two-body interaction is taken to be com- 
pletely random, given by the gaussian distribution of two- 
body matrix elements with V — 0.12 ; this value should 
be compared with the mean energy distance do = 1 be- 
tween the orbitals (single-particle energies) , see details in 
[^,^1^ . The comparison of the thermodynamical temper- 
ature Tfi^ defined by ( |3^ ) with the "canonical" tempera- 
ture (|33) reveals quite strong difference in all the range of 



the rescaled energy x = {E — Efermi)/ {Ec — E 



-^ferrni } 



To 



test our analytical expression for the canonical tempera- 
ture Tcan we have performed direct numerical calculation 
of the temperature according to (||) with the actual spec- 
trum of the two-body random interaction (instead of 
the Gaussian approximation of p{E) ) . Numerical results 
well agree with the analytical expression (|^. 

The knowledge of the equation of state E{T) gives the 
possibility to examine heat capacity of closed systems 
with finite number of interacting particles, 



dE f dA 

dT dE,n 



(37) 



The second term in the above expression is a correction 
which vanishes for highly excited states or in large sys- 
tems, however, it may be important in other situations. 

Following the paper , we can also compare different 
definitions of entropy. Natural definition of the entropy 
in isolated systems can be directly related to the number 
of principal components Npc in exact eigenf unctions., 



(38) 



In such a definition the entropy characterizes the com- 
plexity of a system (note, that for unperturbed "simple" 
states Npc = 1 and Sef = )■ There are several defini- 
tions of Npc one of which is the so-called "entropy local- 
ization length" defined via the information entropy S'inf 
of eigenstates. 



where 



Npc = exp(5inf) (39) 



Sinf{E) = -J2Ek{E) \nFk{E) 



dEkp{Ek)F{Ek - E) \nF{Ek - E) (40) 



Here we have used the i^— function instead of C^^^ 

order to have smooth dependence of the entropy 5'inf on 
the energy E . Another possibility is to find Npc from 
the "inverse participation ratio" , 



(41) 



One more definition is Nz^ ~ max[Ffe(i?)] « F{Ek — E) 
which has been used in 0. The difference between the 
above definitions of Npc depends on a specific shape of 
Fk{E) , however, the values of Npc differ from each other 
by some coefficient which is typically close to 1. 

On the other hand, the estimate for Npc can be ob- 
tained simply from the relation 



(42) 



where V is the spreading width of the function Fk{E) 
and D is the local mean spacing between many-particle 
energy levels. Thus, one can directly relate the number 
of principal components Npc to the density of states p{E) 



Sep = In Npc « In p{E) + In F 



(43) 



One can see that the entropy Se f found from exact eigen- 
states coincides with the thermodynamical entropy Sth if 
the second term in (^) does not depend on the energy. 
As is shown in Appendices 1-3, the spreading width F 
only weakly depends on the energy, in contrast to a very 
strong energy dependence of p{E) . The fact that the 
information entropy S'inf contains the term In p was men- 
tioned for the first time in One should stress that 
the above relations (^, ^3|) are valid if Npc is smaller 
than the size N of many-particle basis, Npc < N/2 . One 
has also to remind that systems under consideration are 
assumed to be in equilibrium, see discussion in Section 
9. 



VIII. INCREASE OF EFFECTIVE 
TEMPERATURE DUE TO STATISTICAL 
EFFECTS OF INTERACTION 

In Section 4 we have shown that in the case of large 
number of particles, the distribution of occupation num- 
bers is of the Fcrmi-Dirac form (^) if the local mean 
field approximation is valid. However, if one uses the 
expression ( p3| ) in order to find the chemical potential /i 
and temperature T , one can obtain inaccurate result. 
To demonstrate this, we have computed the distribution 
of occupation numbers Us for the two-body random in- 
teraction model directly from exact eigenstates of the 
Hamiltonian matrix (|l|) defined in the basis of many- 
particle unperturbed states (see also [^,D). These data 
for the "experimental" values of are shown in Fig. 3 
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by the histogram which is obtained by the average over a 
smah energy window in order to smooth the fluctuations 
(also, an additional averaging over different realizations 
of the random two-body random interaction has been 
done). To compare with the standard Fermi-Dirac dis- 
tribution, we have numerically solved Eqs.(|2|) in order to 
find the temperature and chemical potential. The result- 
ing distribution of the occupation numbers Hs is shown 
in Fig. 3 by circles. One should stress that the value of 
the energy E in ( p3| ) was taken the same as for the ex- 
act eigenstates from which actual distribution of was 
computed, namely, E « i?^*-' . The comparison of the 
actual distribution (histogram) with the theoretical one, 
see (23), reveals big difference for a chosen (quite strong) 
perturbation V = 0.20 . 

This discrepancy is due to the fact that the off-diagonal 
interaction is not taken into account in (|2^). In g it 
was pointed out that one can take into account statisti- 
cal effects of interaction and to correct the total energy 
in (^3|) in the following way. Let us first consider the 
thermodynamical temperature Tj/j defined by (^3|). As 
was mentioned above, see (p^), one can represent the 
width of perturbed density of states p{E) in the form 
= aQ + {AEy where ctq relates to the unperturbed 
density and AE is an increase of the width of the energy 
spectrum due to the interaction. Thus, the temperature 
at a given energy E increases due to the interaction as 
follows 



T = To + AT ; 



{AEy 



Eq — E Ec — E 
resulting in the following relation, 



at/To = {AE)ya'o 



(44) 



(45) 



The explanation of this increase of the temperature 
was given in § and reads as follows. Since the density 
of states rapidly increases with the energy E , the num- 
ber of higher basis states admixed to an eigenstate by 
the interaction is larger than the number of lower basis 
states. An extreme example is the ground eigenstate, 
which contains basis components of higher energies only. 
As a result, the mean energy 



{Ek 



E 

k 



EkF^ 



EkFj;'> po{Ek)dEk (46) 



of the components in an exact eigenstate \i) is higher 
than the eigenvalue i?'*) corresponding to this eigenstate 
(we consider here eigenstates in the lower part of the 
spectrum only). There is another effect which gives the 
increase of (-E/c) j — E^'*-* even if the density of states does 
not depend on the energy. Due to repulsion between the 
energy levels, the eigenvalues move down for this part of 
the spectrum, therefore, the difference between (Ek) and 
increases due to the interaction. This second effect 
shifts the "center" of the function T^*^ = F{Ek - E'-''^) 



. One should stress that all effects leading to the above 
shift of the energy are automatically taken into account 
in the relation (^6|). Thus, one can analytically calculate 
this shift Ae = {Ek)^ - E'-'-^ from the equation (H). For 
this, one needs to know the unperturbed density of states 
and the form of the T— function. The evaluation of the 
shift A E has been done in ^ by assuming some form for 
the F— function which is valid in a wide range of the 
interaction strength V , 



Ae = {Ek 



-Ei^ = ^-^{AEf^^^[AEf 



'-^iE.-E) 



(47) 



Now, in order to obtain the corresponding increase of 
the temperature, one should insert the shifted energy 
E = {Ek)- — + Ae into the equation for the tem- 
perature. 



T = To + AT : 



Er. - E(' 



{AEY 

Er. - Ei' 



(48) 



One can see that Eq.(|48|) is consistent with Eq. (|44|) for 

(A^«a„2. 

Thus, to find correct values for the occupation numbers 
in the Fermi-Dirac distribution, we should substitute the 
increased energy E — E^^^ + Ae = {Ek)i into equations 
( p3| ) for the chemical potential and temperature. The 
resulting shift of the temperature and chemical poten- 
tial leads to the distribution of the occupation numbers 
shown in Fig. 3 by diamonds. As one can see, such a 
correction gives a quite good correspondence to the nu- 
merical data. 

To check the analytical prediction ( p7| ) for the shift 
A^; , we have calculated this shift directly by compar- 
ing the energy T'^'' of exact eigenstates with the energy 
(-Efe),- . The latter has been computed from the exact re- 



lation (T'fc),^ = Efc^fe 



(compare with (p6|)). The 



comparison of these data (circles in Fig. 4) with Eq. 
(straight full line) shows a good agreement, if to neglect 
strong fluctuations around the global dependence. These 
fluctuations are due to fluctuations in the components 
of specific exact eigenstates \i) (note, that the presented 
data correspond to the individual eigenstates, without 
any additional averaging). 

Finally, we would like to note that the described above 
method can be also used to solve the "canonical" prob- 
lem of finding the distribution of the occupation num- 
bers for a system of interacting particles in the thermal 
bath. Indeed, the distribution ns{T) for the "ideal" gas 
is given by the canonical distribution, or more simple, 
by the Fermi-Dirac distribution if the number of parti- 
cles n 1 (in practice, n > 2 is enough). The increase 
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of "kinetic" energy due to the random residual interac- 
tion is given by = {AE)'^/T , see Eq.@. Then, 
one can find the effective temperature which corresponds 
to the same increase of the average "kinetic" energy in 
the ideal gas. It can be done by the differentiation of 
T^^ = dS/dE , where S is the entropy, S = In po + const 
. The result reads as 



In pn , . 
AT = -T——^{AE)^ 
dE^ 



(49) 



For the Gaussian shape of the level density po{E) it co- 
incides with (^). Thus, we can use the Fermi-Dirac dis- 
tribution with the effective temperature T^f f = T + AT 
in order to describe "randomly" interacting particles (for 
AT < T) in the heat bath. 



IX. CRITERIA FOR THE ONSET OF CHAOS, 
EQUILIBRIUM AND THERMALIZATION 

The theory presented above is based on the notion 
of chaos in terms of statistical properties of compound 
eigenstates. Typically, the onset of quantum chaos is as- 
sociated with large number of components in the eigen- 
states. However, as we can see below, this condition is 
not enough for the emergence of equilibrium distributions 
and, in essence, there are different regimes of "chaos". 
Here, we analyze the conditions under which the possi- 
bility of statistical description of isolated systems of in- 
teracting particles can be directly related to statistical 
properties of eigenstate components. 

Let us start with the condition of a large number of 
principal components Npc provided there is an equilib- 
rium distribution of the components of eigenstates. For 
the relatively small interaction, the distribution of the 
components has the Breit-Wigner form. 



F, 



(i) 



1 



PC {Ek-E(^)y + rl^/4 



(50) 



Here, the value Npc = ^ 2D^' defined by the normal- 
ization condition and D = p~^{E) is the mean spacing 
between energy levels, and spreading width is given by 



BW 



27r 



df 



(51) 



where V'^ is the mean squared value of matrix elements 
of the two-body interaction, see (|l]). In the denomina- 
tor df stands for the mean energy spacing of those ba- 
sis components to which a particular basis state \k) can 
"decay" due to direct two-body interaction, see details in 
Appendix 3. Thus, the condition Npc ;» 1 reads 



N„, 



dfD 



> 1 



Vyy^^/djD 



(52) 



However, as was pointed out in [|l8[ (see also 
|l9| , p0|pl]| ) , there is a phenomenon of "localization in the 
Fock space" which means that the states with different 
number of excited particles do not mix with each other. 
This situation occurs when V < df (with some logarith- 
mic corrections). This means that in order to have an 
ergodic distribution for the eigenstate coefficients result- 
ing in the equilibrium distribution of occupation numbers 
for individual eigenstates, one needs both the conditions, 
V > df and (p^). Since df is typically much bigger than 
D , the condition V > df is stronger than (^2|). 

Let us now discuss the properties of eigenstates and 
distribution of occupation numbers in dependence on two 
above parameters, V/df and Npc . Since the value oi Npc 
increases with an increase of V, we first start with a very 
weak interaction for which exact eigenstates have only a 
few relatively large components {Npc ^ 1). In such a case 
the eigenstates are strongly localized in the unperturbed 
basis and, therefore, can be described by conventional 
perturbation theory. This situation is quite typical for 
lowest eigenstates (where the density of states is small) 
even if for higher energies the eigenstates can be consid- 
ered as very "chaotic" {Npc S> 1) . We term this region 
(I) as a region of strong (perturbative) localization. 

The second region (II) is characterized by an "initial 
chaotization" of exact eigenstates and corresponds to a 
relatively large, Npc 3> 1 number of principal compo- 
nents and V < df. The latter condition is essential since 
it results in very strong (non-Gaussian) fluctuations of 
components C^^*' for the fixed energy i?*^'-' of com- 
pound state \i) . Such a type of fluctuations reflects itself 
in a specific character of eigenstates, namely, they turn 
out to be "sparsed" . As a result, the number of prin- 
cipal components can not be estimated as Npc ~ T/D 
, as is typically assumed in the literature. Let us note 
that the energy width F of eigenstates is still close to the 
expression (pl|). 

The above specific properties of compound states can 
be explained by the perturbation theory in the parame- 
ter V/df . Indeed, in the zero order approximation an 
eigenstate \i) coincides with a basis state |fco) for which 
the particles occupy definite orbitals. In the first order 
in V/df the eigenstate \i) is constructed by those basis 
states \ki) which can be obtained from \ko) by moving 
one or two particles (due to the two-body character of 
the interaction). As a result, the coefficients of the state 

\i) can be estimated as cf,*^ = Vk^ ki/{E^'''^ - EkJ ■ If the 
matrix elements Vkg ki are Gaussian variables, for the 

the coefficients C^*'' are also 



fixed spacing — Ek^ 
distributed according to the Gaussian 

The situation is completely different in higher orders. 
For example, in the second order we have C = C'^'^ = 
Efei Vk„k^Vk,kJ{{E^'^ ^Ek,){E'^''^ -Ek,)) and the dis- 
tribution P{C) (for the fixed Is^'^ — EkA ) is close to 
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the Lorentzian, namel y, P(C ) ~ Ccr/C^ for Ccr < C < 1 
where Ccr ~ V^/d'j [g3|,^. Long tails in the distribu- 
tion P{C) are due to the possibility of small values of 
the denominator iJ^'^ — Ek^ in intermediate states. The 
higher orders r of the perturbation theory give stronger 
fluctuations (with additional logarithm in P{C) ) because 
of large number of small denominators. Another type 
of even stronger fluctuations result from very different 
nature of many-body basis states with close energies. 
Namely, the states with nearly the same energy can dif- 
fer by a large number of moved particles in order to ob- 
tain the corresponding states. Therefore, the transition 
into close (in energies) states appear in different orders of 
the perturbation. Therefore, the fluctuations of the com- 
ponents C appear to be abnormal, without obeying the 
standard central limit theorem. In particular, the above 
effect leads to very large fluctuations in the distribution 
of the occupation numbers Tig J QbS cL function of the en- 
ergy E of compound states. Specifically, the fluctuations 
do not decrease as NpJ'^^ , for example, if P(C) ~ 
the fluctuations of (E) do not depend on Np^ at all for 
Npc < I /Ccr (see, for example, [p3[). 

As an example, one can take the two-body random in- 
teraction model (|l|) with n — A particles and m — 11 
orbitals and very weak perturbation V/do ~ 0.02 where 
do is the mean level spacing between single-particle en- 
ergy levels, see Fig. 5. One can see that the "experimen- 
tal" distribution of occupation numbers has nothing to 
do with the Fermi-Dirac distribution (full diamonds), it 
turns out to be even the non-monotonic function of the 
energy of orbitals (see also Q). Note that the av- 
eraging procedure used in Fig. 5 has not wash out the 
fluctuations in . 

With further increase of the interaction, where Npc ^ 
1 and V > df , the region (III) of the equilibriimr for the 
F— distribution (||) emerges. In this region the fluctua- 
tions of the eigenstate components C^*^ are of the Gaus- 
sian form which leads to small fluctuations of the 
occupations numbers in accordance with the central 
limit theorem for the sum (^, Aus/us ^ Npc^'^ <C 1 for 
rig ~ 1 . One should stress that in this region the value 
ol Npc is given by the common estimate, Npc ^ T/D . As 
a result, the distribution of occupation numbers changes 
slightly when changing the energy of a system. Such a sit- 
uation can be naturally associated with the onset of ther- 
mal equilibrium, though the form of the distribution 
can be quite different from the Fermi-Dirac distribution. 
In this case, the F— distribution (||) gives a correct de- 
scription of an actual distribution of occupation numbers 
in isolated quantum systems of interacting particles. One 
can see that the equilibrium distribution for the occupa- 
tion numbers arises for much weaker condition compared 
to that needed for the Fermi-Dirac distribution. Since 
the energy interval df between two-particle-one- hole en- 
ergy levels is small, it is enough to have a relatively weak 
residual interaction V > df in order to have the equi- 
librium distribution (note, that the value of df decreases 



rapidly with the excitation energy, see Appendix 3). 

Next region (IV) is where the canonical distribution (|^) 
occurs; for this case in addition to the equilibrium, one 
needs to have large number of particles, n ^ 1. If, also, 
the condition F <^ ndo is fulfilled, the standard Fermi- 
Dirac distribution is valid with a proper shift of the total 
energy due to the interaction, see Section 8 . Typically, 
this region is associated with the onset of the canonical 
thermalization (see, for example, |p^). 

In practice, the condition (IV) of the canonical ther- 
malization is not easy to satisfy in realistic systems like 
atoms or nuclei since n in the above estimates is, in fact, 
the number of "active" particles (number of particles in 
a valence shell) rather than the total number of parti- 
cles. Thus, the description based on the _F— distribution 
which does not require the canonical thermalization 
condition (IV), is more accurate. 

The above statements are confirmed by the direct nu- 
merical study of the two-body random interaction model 
[^,|| with few particles when changing the interaction 
strength V/do ■ Ifj instead, we increase the number of par- 
ticles keeping the interaction small, V do , the distri- 
bution (|^) tends to the Fermi-Dirac one as it is expected 
for the ideal gas, see 

Finally, we discuss the transition to mesoscopic sys- 
tems. One can show that the result strongly depends on 
the dimensionality d of a system . Let us consider the 
case when the number of particles n is fixed, however, 
the size I of a system increases. Then, the interval be- 
tween single particle energy levels decreases as do ~ ■ 
Since relative interaction between two particles decreases 
likeV ~ I-'' , one can get, V/do ~ r'''-^) . Thus, for 
d = 1 one has V do which means that strong mixing 
(chaos) starts just from the ground state. This is in ac- 
cordance with the absence of a gap in the distribution of 
occupation numbers in ID case (the so-called Luttinger 
liquid). To the contrary, in the 3D case we have ^ ^ do ; 
this means that an admixture of the higher states to the 
ground state can be considered perturbatively, which is 
consistent with the non-zero gap at T = . One can see 
that the transition between regular region (I) and the 
equilibrium (region III) in the 3D case occurs for high 
states when V > df . Recently, the question of this tran- 
sition has been studied in p8| , p0|j2l| . 

X. CONCLUSION 

In this paper we have developed statistical approach 
to isolated finite systems of interacting particles, which 
plays the same role as the canonical approach for sys- 
tems in equilibrium with the thermal bath. It can be 
applied to complex many-body systems like compound 
nuclei, rare-earth actinide atoms, atomic clusters, quan- 
tum dots, etc. The key point of this approach is a new 
kind of the partition function which is defined to the 
shape of compound states (F— function) in the many- 
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particle basis of a system without residual interaction 
(Slater determinants). It allows to calculate mean values 
of different operators as a function of the total energy E 
of a system. As an example, we have calculated the oc- 
cupation numbers ns{E), which may be compared with 
the standard canonical approach giving ns{T) where T 
is the temperature of an open system. In large systems 
(thermodynamical limit) the distribution of occupation 
numbers tends to the canonical distribution with the 
temperature = d [Inp) / dE where p{E) is the energy 
level density. 

Another important area of applications of our ap- 
proach is the calculation of non-diagonal matrix elements 
(transition amplitudes) between the eigenstates of com- 
plex many-body systems. We would like to point out 
that the approach can be also used for solving traditional 
problem of calculations of mean values of operators in 
open systems of interacting particles in the thermostat. 

The suggested approach is entirely based on the sta- 
tistical properties of chaotic compound states which are 
due to the two-body interaction between the particles. 
For relatively strong interaction the number of compo- 
nents of compound states is typically large and these 
components can be treated as random variables, provided 
the two-body interaction matrix elements are "complex" 
enough. 

The essential question is under which conditions the 
above approach is valid in systems with the two-body 
random interaction. Note, the randomness of the ma- 
trix elements itself is not enough for the onset of the 
equilibrium in the system since statistical properties of 
compound eigenstates essentially depend on such param- 
eters as the relative strength of the interaction, excitation 
energy, number of particles and orbitals (single-particle 
states participating in the energy exchange), etc. In par- 
ticular, even if the number of principal components in 
compound states is large, for not strong enough interac- 
tion the statistics of the components can be abnormal, 
leading to huge (non-Gaussian) fluctuations in the struc- 
ture of the eigenstates. In such a case, there is no equilib- 
rium in the system and standard statistical description 
is not valid. In this regime, our numerical data show 
that the distribution of the occupation numbers strongly 
fluctuates when slightly changing the total energy of the 
system. Therefore, the transition to the statistical equi- 
librium is far from being trivial in the systems with finite 
number of interacting particles. However, for larger in- 
teraction the fluctuations of the eigenstate component 
becomes normal (Gaussian) and the equilibrium occurs. 
In this situation the fluctuations of the occupation num- 
bers ns{E) are relatively small and the i^— function gives 
a correct result for the isolated finite systems even for 
small number of interacting particles (provided the num- 
ber of principal components in the eigenstates is large) . 

The advantage of our approach, in comparison with 
direct calculations of complex quantum systems, is that 
we do not need to diagonalize huge Hamiltonian matri- 
ces in order to perform calculations of observables for 



excited states. Indeed, for the full statistical description 
of such systems, one needs to know the average shape of 
the compound states (rather than the eigenstates them- 
selves) and the unperturbed energy spectrum. There- 
fore, the problem of analytical description of the shape 
of eigenstates is the central point in the technical imple- 
mentation of the approach. One should note that the 
average shape of eigenstates is the same as that of the 
local spectral density of states (LDOS), if the interaction 
is not extremely strong. 

For small interaction the shape of the chaotic eigen- 
states is known to be well described by the Breit-Wigner 
form. However, in practice, this region is small if the 
number of particles is not very large. With an increase of 
the interaction strength, the average shape of the eigen- 
states (F— function) changes from the Breit-Wigner one 
to that close to the Gaussian with the exponential tails. 
As we have found (Appendix 1), the correct description 
of the shape requires two essential parameters. The first 
one is, in fact, the half-width of the F— fimction which is 
close to the Breit-Wigner half- width and for weak inter- 
action is given by the Fermi Golden Rule. Another pa- 
rameter is defined by the root-mean-square width of the 
F— function, or, the same, by an effective band-width of 
the Hamiltonian matrix in the energy representation. 

In this paper we suggest the phenomenological expres- 
sion for the F— function which is valid in a large region of 
the interaction strength and other parameters. This ex- 
pression allows for the analytical and numerical calcula- 
tions of different mean values and transition amplitudes. 
As one example, we have derived analytical expression for 
the distribution of the occupation numbers ns{E) in iso- 
lated systems of n interacting Fermi-particles distributed 
over m orbitals (we have used the simple Gaussian ap- 
proximation of the F— function, see Section 5). 

By making use of the F— function, we have also stud- 
ied the validity of the standard Fermi-Dirac distribution 
for the description of finite systems of interacting Fermi- 
particles. As was found, the Fermi-Dirac distribution 
can provide reasonable approximation for both isolated 
and open (in the thermal bath) systems. However, the 
parameters of the Fermi-Dirac distribution have to be re- 
defined by taking into account the increase of the effective 
temperature which is due to the effects of the interaction. 
We calculated this increase of the temperature analyti- 
cally and compared with the numerical experiments for 
the two-body random interaction model; the data show 
very good agreement. 

One should stress that our approach gives more ac- 
curate result for rig and has much wider region of the 
applicability, compared to the Fermi-Dirac distribution. 
Specifically, it is valid even in the region where the Fermi- 
Dirac distribution fails, for example, due to the small 
number of particles. 

To conclude with, we would like to point out that sim- 
ilar approach may be also used for classical chaotic sys- 
tems. Indeed, let us consider the system described by 
the Hamiltonian H = Hq + V where Hq is the unper- 
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turbed ("simple") Hamiltonian and V stands for "com- 
plex" interaction between particles. Assume that we 
know the distribution for some variable in the system 
described by Hq , for example, the statistical average 
for the energy distribution of a single particle n(e, Eq) 
where e is the energy of the particle and Eq is the 
total energy of the system. Then, one can calculate 
the effect of the interaction V by averaging of n(e, Eq) 
over the unperturbed energy Eq using the -F— function, 
n{e,E) = Jn{e,EQ)F{E,EQ)dEQ. Here F{E,Eq) gives 
the probability of different values of Eq for a given value 
of E. As was indicated in and checked numerically in 
1^,^, this classical F— distribution coincides with the 
shape of quantum eigenstates in the semiclassical region 
which turns out to be very wide. Thus, the knowledge 
of the shape of quantum eigenstates can be used for the 
classical calculations and vice versa. 
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Appendix 1. Structure of chaotic eigenstates 
and spreading function 

For practical implementation of the above approach 
one needs to know the average shape of compound eigen- 
states (the F— function). One should stress that there 
is no simple analytical expression valid in a large range 
on the interaction strength V . For example, the pop- 
ular Breit-Wigner expression (Lorentzian) is not good 
for obvious reason: it has infinite second moment. The 
question of an appropriate description of chaotic eigen- 
states in realistic many-body systems has been studied 
in detail for Ce atom In particular, it was found 
that good correspondence to the numerical data is given 
by two phenomenological expressions. The first one is 
F{x) ~ exp(-Vl -I- 4x2) wiiere x ={Ek- E)/T with Ek 
as the energy of a basis state |A;) and T as the effective 
width of the distribution. This expression is close to the 
Gaussian at the central part and is exponential in the 
tails (see similar conclusions in | [T^ where nuclear shell 
model was studied). 

Another expression which is more convenient for the 
analytical study is the so-called "squared" Lorentzian , 



Eh — H, 



F{Ek - E) 



1 



{Ek - EY 



£2 

4 



E = F(*) + Ai' 



(53) 



Here A^*'' <C F is some small shift (see below) which 
in the zero approximation can be neglected, and E^ is 
defined by 



kk 



s>p 



(54) 



Since the "resonant" dependence E^ — E'^^^ of the 

(i) 

spreading function F^ for not extremely strong inter- 
action is symmetric in indexes i and k , the value of 



T = 2\l (AF) can be expressed in terms of the second 
moment of F using the following exact relation for the 
basis components, see Appendix 2, 

2 / , ,\ 2 



(AF)^ 
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p^k 
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kp 



(55) 



with Hkp standing for non-diagonal Hamiltonian matrix 
elements defined by the residual interaction V . This al- 
lows us to find the second moment of the spreading func- 
tion F(F). For example, in the case of n particles dis- 
tributed over m orbitals we have [pi 



- ^(AF)^ = — n(n 



l){m — n)(3 + TO — n) 



(56) 



Here — \Vst^pq\ is the mean squared value of non- 
diagonal matrix elements of the two-body residual inter- 
action. 

Our detailed study of the two-body random interaction 
model has revealed that the shape of the eigen- 

functions. (as well as the local spectral density of states) 
strongly depends on the relative strength of the inter- 
action. Namely, with an increase of the interaction V , 
the shape of the F— function changes its form from the 
Breit-Wigner one to the nearly Gaussian (see also |13|). 
It was found that for small residual interaction the shape 
of eigenstates has more complicated form, compared to 
(^), and should be characterized by two different widths. 
Indeed, the half-width of the F— distribution is given by 
the Fermi golden rule Tbw — 2'Kv'^/df where v is the ma- 
trix element of the residual interaction coupling a partic- 
ular basis component with other basis states |/) directly 
coupled by the two-body interaction, and df is the en- 
ergy spacing between these basis states (see details in 
Appendix 3). On the other hand, there is the relation 



( p6| ) which defines another width F = 2\J (AF)^ via the 
second moment. One should stress that these two widths 
are parametrically different in the interaction, F^^y 
, and T ^ V . There is also the "non-resonant" energy 
dependence F cx p^^ (slow variation of the F— function 
due to the change of the density of states p{E)) which 
should be taken into account. This dependence follows 
from the estimate F^ax ~ ^ Tp. 

The above arguments allow us to find more universal 
expression for the spreading function F when Tbw < T 



F 



(i) 



(po(i?;c)p(FW)) 



-1/2 



{Ek - E)' 
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(57) 
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Here we take into account the shift of the maximum of the 
i^— function by the relation E = ii^^'^ + A^*^ . The two pa- 
rameters, Fi and r2 are directly related to the above two 
widths, Ti =Tbw and T2 = T^/Ti. The value of T2 is 

found from the relation J2^ Pjc^ (E^'^ - Ekf = JaW = 
(see ( p^ ) and Appendix 2) by integration of ( p7| ) in 
the approximation p = const. Here and below we assume 
that the F— function is normalized, Z = 1. In the expres- 
sion (56), po{Ek) is the density of basis (unperturbed) 
is the density of compound states. 



Po{Ek} 

states and piE^"^- 
We assume they are smooth functions of the energy and 
this energy dependence is slow in comparison with the 
"resonant" energy dependence on the scale F. Symmet- 
ric dependence on po and p has been chosen in order to 
keep symmetry in indexes k and i in the function. 

For very small V we have Fi ^ F2 (also, A^*' ^ 
Fi, see below) .therefore, in the central part the 
i^— distribution ( pT\ ) has the Breit-Wigner shape with the 
width Tbw ■ Concerning the meaning of F2 , it is the 
effective energy band width of the Hamiltonian matrix H 
, which is due to the two-body nature of the interaction. 
Indeed, the expression for F2 is given by the estimate 



BW 



Of 
27r 



n{n — l)(m — n)(m — n + 3) ~ dQ{m — n) 

(58) 



which is independent of the interaction strength V . Here 
we have used expression (^) and the estimate of the av- 
erage value df = do/Mf for high excited states {Mf here 
is the normalized density of those basis states which are 
directly connected to the chosen state, see details in Ap- 
pendix 3). On the other hand, the typical band width 
Ah = 2bdo of the two-body interaction Hamiltonian ma- 
trix is about four times. Ah ~ 4(m — n)dQ , of the en- 
ergy needed to transfer the particle from the Fermi-level 
ep = ndo to the highest available orbital = mdo ■ 
Therefore, the estimate for F2 reads as F2 ~ A///4 . 

Now, we can easily explain the form of the i^— function 
(^) using the perturbation theory in the interaction 
V . First, let us consider the energy interval Tbw < 
Ek — E^^^ I < Ah ■ Within this interval, the basis state 
k) can be coupled to the principal components of the 



( V 



stateji) in the first order of V 

. This quadratic decay agrees with the Breit-Wigner 
shape of the i^— function. Outside the energy band Ah 
, for Ah < \Ek — -E^*-*] < 2A^f , the basis state \k) can 
be coupled to the principal components of the state |i) 
in the second order of V , resulting in the dependence 



This corresponds to the tails of 



the squared Lorcntzian shape (|57D. Therefore, our ex- 
pression ( |57| ) seems to be good in a large energy interval 
and the second moment is finite which is important for 
applications. And finally, longer tails are described by 



higher orders v = - — ^ of the perturbation theory, 



(i) 



Ek - E('') 
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(59) 



This explains the exponential tails of the _F— function, 
see details in 

Numerical calculations demonstrate that 

for stronger interaction V , the width of the spreading 
function F rapidly becomes linear in V (instead of the 
quadratic dependence in Tbw ) and it is better to use 
(p3|). One can write the extrapolation expression both 
for small and large values of V (see also |lq]): 



Fi 



BWl 



Tbw + F 



(60) 



As a result, for small V we have Tbw ^ F and Fi = 
Tbw ~ and for larger values Fi « F2 w F 1/ > 
A^*-* . The critical value for this transition is given by 
the relation Tbw = F = F2 ^ Ah and reads as Vcr ~ 
dfn{m — n)/(27r) , see also (p6|). The estimate of the 
average value of df = d^/Mf (see Appendix 3) far from 
the ground state gives Vcr ^ do/n . As was discussed in 
Section 9, the equilibrium distribution occurs for V > df 
, this results in a quite unexpected conclusion. Namely, 
the validity of the standard Breit-Wigner shape turns out 
to be very strongly limited since the region 1 < V/df <C 
Vcr/df ~ n(m — n)/(27r) is practically absent for small 
number of particles n and orbitals m . 

The shift A^'' in ( ^ ) stands due to the level repul- 
sion which forces eigenvalues i?*-*' in the lower half of 
the spectrum to move down. The mean field energies 
Ek = Hkk do not include the non-diagonal interaction 
which leads to the repulsion, therefore, the "center" of 
the F— hmction is shifted by the value A^*^ — Ha — i?^') 
where Ha is the diagonal matrix element of the Hamilto- 
nian matrix. This shift A^*^ can be estimated from gen- 
eral arguments. Indeed, the shape of the density of states 
is the same for both interacting and non-interacting par- 
ticles , with the same position Ec of the centers of 
Pq{E) and p{E) (due to the conservation of the trace of 
the Hamiltonian H ). However, the variances of po and 
p are different. This means that one can use the scaling 
relation {D KD) for the energy intervals D and find 
the scaling coefficient K from the relation between the 
variances, cr^ = o-q^ + (AE)^ = K'^ctq where (AE)^ is 
defined by (p6|). Since the center Ec for the energy level 
density p{E)does not shift, one can obtain the following 
shift of the levels: 




(61) 
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The value of (AE')^ is typically much less than ctq^ , 
therefore, one can get 



2(ao)2 



(62) 



Another way to obtain the shift A^'-* is related to the 
exact relation for the first moment of F , see Appendix 
2, 

Ek =J2E^'^Ek^ ~ J fI;^ piE^'y)E^'UE^'^ (63) 



The substitution of the expression (^ 
in the following value of the shift ||] , 



into (p3) results 



(64) 



According to ||6|,[7| (see also ||11|) the shape of density of 
states for m ^ n ^ 1 is close to the Gaussian both for 
non-interacting and interacting particles with Ec and ctq 
as the center and the variance of the energy distribution 
Pq{E) (respectively, tr^ for p{E)). In this case the relation 
(U) gives the same estimate (^) for the shift A^*"*. 

The fact that the two different derivations of the shift 
A[*^ lead to the same result is far from being trivial since 

the assumptions for the two derivations of aS*^ are dif- 
ferent. Indeed, the second derivation of (^ is based 
on the specific dependence of the eigenstate shape on 
the densities po and p, unlike the general derivation of 
(|6l|). One should stress that the specific form of the 
"resonant" energy dependence of the F— function (the 
denominator in (|5^) defined by the squared Lorentzian, 
Gaussian, etc.) is not important for Eq. ( |64[ ) provided 

(Aii')^ is fixed. The only assumption in the above deriva- 
tion is the possibility to expand the density p{E) near 
the maximum of the i^— function. In fact, above we have 
demonstrated that the non-resonant "distortion" factor 
i ={pQ{Ek)p{E^''^)y^'^ in d^) is necessary. 

Thus, the phenomenological expression for the shape of 
the F— function ( p7| ) is self-consistent. Note, that one can 
use other expressions for the F— function (see, for exam- 
ple, [|3|), however, it should contain both the resonance 
term depending on Ek — E'^^^ and the density distortion 
factor ^ . 

Appendix 2. Moments of the F— function and 
energy spectrum 

Here we calculate the first and the second moment 
of the function F'jf' over the perturbed spectrum F^*) 

. Note, that the dependence of F^'^ on the energy F^*) is 
known as the local spectral density of states (LDOS), or 
the "strength function" . On definition, 



\H\j) {j\k)^Hkk^Ek 



(65) 



where the relation (i| H \j) = 5ij {i\ H \i) is used for the 
exact eigenstates. The variance can be obtained using 
the matrix elements of 



Ek - F(^ 



(66) 



For example, in the two-body random interaction model 
with n particles distributed over m orbitals the sum in 
( |66| ) can be evaluated exactly (see (p6|)). 

Now, we calculate the first moment and variance of the 
energy spectrum. The trace conservation of H gives the 
first moment, 

F =-^H -(-^H V 

c kk I ^ kk I 

k \ k J 



The conservation of TrH^ gives 

^ = ^E(^i^'i^) = 4Ei^^: 



which results in the relation 
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where ctq is the variance of the unperturbed spectrum 
and we have used Eq.(66). 

Appendix 3. Calculations of spreading widths 

To start with, we should stress that there are different 
definitions of the spreading widths. One of the natural 

definitions is Tfc = 2J{AE)l where (AF)^, is the vari- 
ance of the distribution of the components F^*'' (see (|66|)). 
In Ref. it was shown that in the model of random two- 
body interaction this quantity is constant, Ffc = F , see 
(p6|), i.e., it does not depend on a particular basis state 
(therefore, on excitation energy, number of excited par- 
ticles n* corresponding to this state, etc.). Note, it has 
a linear dependence on the interaction strength V. 

Unlike the latter, the commonly used definition of the 
spreading width is related to the Breit-Wigner distribu- 
tion and is defined as its half-width Tbw ■ However, this 
definition is reasonable only for relatively small interac- 
tion, when the form of the F— function is indeed close to 
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the Breit-Wigner form. In this case, the spreading width 
is given by the Fermi Golden Rule, 



BW 



2-K- 



if 



(67) 



where is the mean spacing between corresponding ba- 
sis states I/) and 14/ is the matrix element of the inter- 
action between the basis states |fc) and |/) Q. As one 

can see, the spreading width is proportional to V'^ 

and differs from F . Note, that the second moment of 
the Breit-Wigner shape diverges, however, actual form 
has always cut-off in the tails which is reflected by the 
finite value of F . 

If the interaction V is not small, the form of the 
i^— function significantly differs from the Breit-Wigner 
shape (see Appendix 1). The critical value of V for this 

transition can be estimated from the condition F 



{k) 

BW 



Tk ■ 

Contrary to the spreading width F , the half-width of 
the Breit-Wigner shape depends on the basis state |fc) . 
Let us start with the estimate of the mean value of the 



spreading width F^^ . According to the definition (67), 
one needs to calculate the density of squared transition 

matrix elements V^^/df ~ (X)/ ^fc/) where the sum 
is taken over transitions from a given basis state \k) to 
other basis states |/) in the energy interval A . One 
should stress that the number of "allowed" transitions 
(due to the two-body interaction) is much less than the 
total number of basis states in this interval A . Let 
us assume, for simplicity, that the spacings do between 
single-particle energy levels is constant and consider a 
pair of particles which occupy the orbitals s and q . If we 
move one particle to higher orbital and another particle 
to lower orbital by the same energy interval, the total 
energy of particles does not change. The total number 
of such moves is M « (min(TO — s,q) + min(TO — q, s))/2 
where m 3> 1 is the number of orbitals. By averaging 
over all values s and q , one can obtain M m/S . 
The Pauli principle reduces the number of available or- 
bitals, therefore, in the system with n particles we have 
M « (to — n)/3 . The number of possible pairs is given 
by n{n — l)/2 , thus, the total number of basis states 
I/) which have the same energy and connected with the 
chosen basis state \k) is defined by 



Mf K. {m~ n)n{n - l)/6 



(68) 



Other basis states are separated by, at least, the energy 
distance dp . As a result, we obtain 



BW 



27rMf — Ri -(m-n)n(n-l) (69) 

do 3 rfo 



It is interesting to note that the simple estimate in- 



volving total strength of transitions (AE) , see Eq.(56), 
divided by the band- width A ^ do{m—n) , gives close re- 
sult. The above estimate (p9) can be used when studying 



the shape of the F— function in the regime of weak inter- 
action. Note, the value of Mf defines the mean spacing 
between the "allowed" final states is 



df = do/Mf 



(70) 



One should stress that the above estimate of Mf has 
been obtained for the spreading width averaged over all 
basis states. Near the ground state actual value of Mf 
is much smaller due to limitation of the available phase 
space. Also, the spreading width depends on the number 
of excited particles in the basis state. For example, one 
can calculate the spreading width Tg^, of the basis states 
with one excited particle only. 



Mf 



F^ 

^ BW 



do 



s<2n (71) 



M/ = n(2s-3n); F^ - 3")^o' . ^ > 2n 

do 

(72) 

Here s is the position of a particle corresponding to the 
energy w s do ■ 
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FIGURE CAPTIONS 

Fig.l. Analytical description of the occupation num- 
bers. The data are given for the two-body random in- 
teraction model (|l|) of n = 4 Fermi-particles distributed 
over m = 11 orbitals with V = 0.20 and do = 1 in the 
definition of single-particle energies, = dQ{s + i), see 
|^,||j|]. The histogram is obtained according to @j by 
the averaging over eigenstates with energies taken from 
small energy window centered at E — 17.33 and over 20 
Hamiltonian matrices (|l|) with different realization of the 
random interaction. Stars represent the analytical ex- 



pression (28) with aos found from single-particle energy 
spectrum. Diamonds correspond to the Fermi-Dirac dis- 
tribution with thermodynamical temperature (|3^) and 
chemical potential found from the standard condition for 
the total number of particles, n = Us ■ 

Fig. 2. Different temperatures versus the rescaled en- 
ergy {E - Efernii)/{Ec - Efernii) for the two-body 
random interaction model with n = 4 Fermi- particles 
and m = 11 orbitals. Triangles stand for the thermo- 
dynamical temperature Tth defined by (33) and should 
be compared to the canonical temperature Tcan (circles), 
see (Q). The width a of the perturbed density of states 
is defined by the residual interaction V — 0.12 accord- 
ing to (^8|) and ( p6| ) with ctq found numerically from the 
unperturbed many-particle energy spectrum (the mean 
level spacing do between single-particle levels is set to 



do = l). 

Fig. 3. Fermi-Dirac distribution for strongly interacting 
particles. The data are given for the two-body random 
interaction model (|^) with the parameters of Fig. 1 (the 
rescaled energy is x = 0.55). Circles stand for the Fermi- 
Dirac distribution with the total energy E corresponding 
to the energy of eigenstates, see (^3|). Diamonds corre- 
spond to the shifted energy according to the expression 

Fig. 4. Shift of the total energy for the corrected Fermi- 
Dirac distribution. The data are given for the model (|l|) 
with n = 4, m = 11, do = 1,V = 0.12, see explanation 
in the text. The straight line is the analytical expres- 
sion (|4^); the dotted line (circles) present direct compu- 
tation of the shift based on the diagonalization of the 
Hamiltonian (^) with the following computation of the 
< Ek >i ■ On the horizontal axes the rescaled energy 

X(') = (£;(') - Ef,rmi)/{Ec - Eferm^) is plotted. 

Fig. 5. Distribution of the occupation numbers for 
small interaction. The histogram is obtained in the same 
way as in Figs. 1-2, for the very weak interaction V = 0.02 
which correspond to the region of (II) of the "initial chao- 
tization", see Section 10. The total energy (center of the 
small energy window) is -E'*' = 17.33. Diamonds corre- 
spond to the theoretical expression (^ which is not valid 
in this region due to absence of equilibrium. Stars are ob- 
tained by direct numerical computation of Us according 
to the F-distribution (||) with the F— function taken in 
the form (57), see Appendix 1 and ||^. The latter values 
are closer to the "experimental" ones since we performed 
summation over real unperturbed spectrum (instead of 
the integration with the Gaussian approximation for po 
used to derive Eq. (|2|)). 
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